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Abstract. - The change of the structure of concentrated colloidal suspensions upon addition 
of non-adsorbing polymer is studied within a two-component, Ornstein-Zernicke based liquid 
state approach. The polymers' conformational degrees of freedom are considered and excluded 
volume is enforced at the segment level. The polymer correlation hole, depletion layer, and 
excess chemical potentials are described in agreement with polymer physics theory in contrast 
to models treating the macromolecules as effective spheres. Known depletion attraction effects 
are recovered for low particle density, while at higher densities novel many-body effects emerge 
which become dominant for large polymers. 



Colloid-polymer mixtures (CPM) play an important role among dispersion systems for 
quite different reasons. On the one hand their technological use has long been realized and 
the specific colloid-colloid interaction caused by free polymer, the depletion attraction 
is exploited to induce flocculation or phase separation in dispersions. On the other hand, this 
system allows one to experimentally address the fundamental question about the requirements 
on the pair potential for a liquid phase to exist 0,0]. For even more complex systems CPM 
serves as a model system in order to address protein crystallization and other phenomena 
involving spherical nanoparticles. 

Recent field theoretic considerations of a few colloidal particles in dilute polymer solu- 
tions |(| , and scaling arguments for semidilute solutions fj| , have clarified the polymer struc- 
ture close to particles and provided a fundamental understanding of the origins of the depletion 
attraction. The phase diagram of sterically stabilized hard-sphere like colloidal particles in 
polymer solutions close to their O-point has been mapped out in detail jj| , and a rather suc- 
cessful mean-field like theory for it exists M . The thermodynamic approach of Lekkerkerker 
et al. maps the CPM onto a binary hard sphere mixture with non -additive radii, allowing the 
effective polymer spheres (EPS) to overlap but excluding them from the colloids. Simulation 
studies for a closely related model support some of the theoretical predictions || . 

So far little, however, is known about the structure of the CPM, nor how depletion phe- 
nomena change when particles are small compared to polymers. Whereas neutron scattering 
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experiments [ fiO| on rather low density colloids could be fit into an effective pair potential 
description, more recent light scattering experiments on more concentrated systems |]TT| have 
evaded theoretical explanations with EPS models and the assumption of a colloid struc- 
ture unperturbed by polymer fails P,pj]]. Removing the assumption of an unaffected local 
colloid structure forces one to develop a new approach to CPM which explicitly addresses 
local structure. Moreover, as the packing of polymers into the void space between particles is 
of intrinsic interest on its own, a microscopic two-component approach appears desirable and 
shall be presented in this letter. This macromolecular approach further provides the unique 
possibility to explore the full range of colloid-polymer size ratios. 

We consider a two-component system consisting of hard spheres of diameter a c at packing 
fraction cf) c = ^g c a^, and polymers modeled as chains of segments (excluded volume diameter 
cr p ) characterized by a (Pade-approximated) Gaussian intramolecular structure factor oj(q) 
|13| , where q is the wave vector. In obvious two-dimensional matrix notation, with a diagonal 
matrix of intramolecular structure factors (ujh = uj{q), and LO22 = 1), diagonal matrix of 
site number densities, Qij, and a matrix of direct, dij(q), and intermolecular, hij(q), site- 
site correlation functions, the Ornstein-Zernicke-like equations for the total structure factors, 
S l3 (q), are 

S-^q) - [gu>(q) + Qh(q) Q \- 1 = V 1 - %) • (1) 

Together with the constraints of (additive) local steric exclusion, gtj(r < |(cr, + <jjY) = 0, 
where the gij(r) = 1 + hij(r) are the intermolecular pair correlation functions, eq. (|lj) may 
be viewed as a definition of the effective interactions Cij(r). For the pure colloid component 
we adopt the well established Percus-Yevick (PY) approximation, c cc (r > a c ) = 0. For the 
pure athermal polymer component, the polymer reference interaction site model (PRISM) 
approach is employed as it has proved versatile and successful for polymers Jll|[l4]]. PRISM 
considers the segment averaged polymer structure and, in the most simple version, adopts 
a PY-like closure, c pp (r > a p ) = ]l3 , 14 [ , also derivable from a Gaussian field theoretic 



perspective |lf_ 

To close the integral equations, a further approximation for the colloid-polymer direct 
correlation function is required. We propose a novel generalization of the PY closure motivated 
by the known physical behavior of polymer packing near a spherical particle [|| : 

On the segmental level, again the central idea of a short-ranged steric interaction is used. But, 
an effective entropic repulsion between segments and particle is present over a length scale A 
beyond contact due to chain connectivity constraints. On physical grounds A cannot be larger 
than the smaller of the two sizes (colloid or polymer). In general, the nonlocality length A will 
depend on the densities and relative sizes, and additional information is required to determine 
it. We use the argument of thermodynamic consistency and calculate the excess free energy 
of insertion per segment of a polymer coil into the pure colloidal hard sphere solution via two 



independent routes |14|. From the compressibility route of concentrating the colloid solution 



around the polymer {(3 = l/fe^T): 



Be 



PS^\ ep=0 = dg' c c cp (q = 0, q' c ) , (3) 

Jo 



and from the charging route which considers the free energy required to grow the particles 
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Fig. 1 Fig. 2 

Fig. 1 - Excess chemical potential N/35f/!p\ g =o for inserting one polymer into a hard sphere fluid 
versus polymer size for three colloid packing fractions. The present results, where error bars mark the 
differences from eqs. (^[^), are compared to the equivalent PY results for inserting a sphere of radius 
equal to R g (dashed lines) as used in the phantom sphere approach of Ref. |]]. The inset compares 
the excess free energy cost (solid line from Eq. ^ where fig = jC p cr^/P) for adding one colloidal 
sphere to a dilute polymer solution as a function of relative polymer size, £ c = R g /(s/2cr c ), with the 
field theoretic result, dashed line Q. 

Fig. 2 - Spinodal curves for various size ratios. The polymer size independent asymptote (f) c — ► 1/22 
for tp p — > oo is marked by a vertical bar. The limits for small or large polymer sizes are included as 
thin solid or dot-dashed line, respectively. 



from points to their actual sizes: 

P5^\ ep=0 = Jdaa p ^a c fg^C-^) + Shrfifo ^ C (^c) 2 ^|^ | g ^o -(4) 

Equating the excess free energies from both routes determines A, and, as we first focus on dilute 
polymer solutions, determining it at g p = appears justified. As a technical approximation we 
analytically determine the leading limiting behaviors and match the following Pade-form to A: 
A -1 = + Jur^' wnere — (\/5— l)/4. Here £ is the (collective) polymer correlation 
length and determines the width of the depletion layer in the large particle limit fj, 13|. For 



accessible colloid densities and size ratios, the different routes, eqs. to the insertion free 

energies at vanishing polymer concentration then agree with relative errors smaller than 15 
%. Representative results are shown in fig. |]. Using the results for £ from eqs. (|l|) to ([|) 
at finite polymer concentrations provides a reasonable approximation to the exact A at all 
densities. The chemical potential to add single colloidal spheres to a polymer solution from 
the compressibility and free energy routes then agrees within an error of a factor 3 at most. 

The aspect of a large length scale separation, a c /a p w 10 3 , thwarts numerical studies 
of many realistic systems, but can be exploited in a scaling limit of shrinking the polymer 
segment size, a p , to zero; for a rigorous derivation of this field theory inspired "thread-limit" of 
PRISM for Gaussian chains see [^6|. Then the system is characterized by the colloid packing 
fraction, (f> c , the reduced polymer concentration, ip p — 2ttc p ^, where c p is the polymer- 
molecule number density, and the polymer size £ c = R g / y/2, where R g is the radius of gyration. 
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Note, that we now use the colloid diameter as unit of length; er c = 1. Flexible Gaussian 
polymers are obtained where the statistical segment size \[%o p follows from eqs. (||j4|) for 
q p — > and g c — > 0. The effective polymer concentration, ip p oc c p /c*, is measured relative 
to the dilute-semidilute crossover concentration, c*, where the polymer coils start to overlap. 
This first result from the present macromolecular approach justifies EPS models, where rj p = 
^-CpRg w </? p /0.53 is used from the outset. For <p c = 0, the polymer correlation length is 
£ = £ c /(l + 2<y9 p ) |u|. The connection of the thread limit of PRISM to Gaussian field theory 
and scaling approaches has been discussed Jl3|.p^|. 

The non-linear integral equations can now be solved with the Wiener-Hopf factoriza- 
tion technique; details will be presented elsewhere. The lines of spinodal instabilities, where 
the (partial) compressibilities diverge indicating fluid-fluid phase separation, can be found 
accurately without the need of extrapolations of numerical integration schemes. In fig. ||, 
one observes the trend that the required dimensionless polymer density increases as R g grows 
and/or particle diameter decreases. This prediction |l7|] is in qualitative agreement with many 
experiments on polymer- colloid, protein or micelle suspensions but in disagreement 

with EPS models which do not account for polymer-polymer repulsive interactions, many 
body depletion effects, nor particle penetration of polymer coils, all of which tend to reduce 
the tendency for phase separation. 

The polymer segment profile close to a colloidal particle is correctly predicted to be of a 
parabolic form: g C p(r)=-^{r — a c /2) 2 + . . ., where A — > | + (A+£)/cr c for <fi c — > 0. For nonzero 
A, and simplifying to small polymers, £ c <C cr c , the number of polymer segments in contact 
with the colloid particle is of order unity for a single polymer and becomes independent of the 
polymer degree of polymerization in the semidilute regime, g p g C p(r ~ <Tc ^ crp ) oc c p (l+2co p ) 2 — > 
(g p a p ) 3 . This behavior is in agreement with scaling Q| considerations. Quantitatively we find 
good agreement as ip p — > 0: A — ► 1/2 or A — ► |r for £ c — * or £ c — * oo, respectively, 
compared to A — > 1/2 or A — > 4Ai| SL from field theory for a single Gaussian chain ||. For the 
PY closure, where A = 0, this behavior is violated and thus the polymer induced depletion 
attraction is significantly underestimated . The free energy of insertion of a sphere into a 
polymer solution can be obtained from the free energy route (analog of eq. (Q)) as: 



/^L=o = + m 5 ^) + f (1)2) ■ (5) 

6£ z 6Ai <T C M <Jc 



For dilute polymer solutions, eq. (|^), is compared with field theoretic calculations |6) in the 
inset of fig. [I], and excellent agreement is found. Whereas the colloid excess chemical potential 
measures the number of excluded polymer molecules for small coils, for large polymers, £ 3> <r c , 
only segments along a strand of length proportional to the colloid size need to be rearranged, 
and the result becomes intrinsic (independent of polymer size), /35/j, c \ gc =o oc g p a c a p . For 
identical physical reasons as in eq. (|^), the chemical potential for inserting a large polymer into 
a sphere fluid becomes intrinsic, f3S^, p \ gp= o oc R X J U /N oc N° with v = i, as the small spheres 
need to accommodate (D(Ry v ) polymer segments. This asymptotic behavior is apparent in 



fig. |l| and qualitatively agrees with RISM-based theory for a single delocalized electron in 
a hard sphere fluid (§0|. EPS models, which predict (35^ p \ ep=0 oc R 3 g /N oc N 1 / 2 @, already 
begin to deviate from PRISM for £ c < <t c , as seen in Fig |[ This iV-dependent overestimate 
of the polymer insertion chemical potential has direct consequences on the predicted trends 
for phase separation as a function of Cc/c c , which are opposite to our results in fig. |^. 

Figure [?] shows the pair correlation functions for different parameter sets with increasing 
colloid concentrations. For rather short polymers, fig. ||(a), the polymer induced depletion 
attraction leads to a strong increase of the colloid contact value, g cc (r — cr c +), relative to the 
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hard sphere value. The polymer depletion layer is visible in g C p(f) for r — h < £. Outside of 
the depletion region oscillatory correlations are found at higher <p c reflecting the imprinting 
of colloid order on segment-particle packing. Clustering results in the interchain polymer 
segment correlations, g pp (r), and fills the correlation hole, which in pure dense polymer flu- 
ids results from effective intermolecular repulsions extending out to £ c fTs}] . As the polymer 
correlations may be considered to decay for r < a c , the description of the colloid structure 
when <7 C S> £ c appears amenable to traditional one-component effective pair potential ap- 
proaches pi pi, which rely on the assumption of negligible long-ranged polymer correlations. 




Fig. 3 - Pair correlation functions for two polymer sizes £ c . In (a), £ c = 0.03 and reduced polymer 
concentration ip p — 0.061, for increasing colloid packing fractions, (f> c = 0.05, 0.2 and 0.4; <? cc (r) 
increases to the contact values of 5.80,8.11 and 14.1. In (b), £ c = 0.202 and identical (f> c are shown 
at ifip — 0.086. For both polymer sizes the small-q limit of the colloid structure factor agrees at 
(f> c — 0.40, where S cc {q = 0)/g c = 0.80. The thin solid lines in both panels labeled HS give the 
PY-result for hard spheres at <j> c = 0.40. 

For larger polymers the depletion layer in fig. ||(b) widens as expected. Moreover, increas- 
ing colloid volume fraction leads to much longer-ranged polymer fluctuations as strikingly 
apparent for g pp at <p c = 0.40. The clustering enforced by the void structure among the col- 
loidal spheres overwhelms the polymer repulsion responsible for the correlation hole. As the 
polymer correlations extend beyond a number of particle diameters, effective pair potentials 
cannot be used. Three- and higher body contributions to the effective colloid interactions be- 
come important which are mediated by the increasingly correlated polymer structure; this is 
in qualitative agreement with simulation studies of a simplified model using interpenetrating 
polymers (61, which should apply to our systems for low polymer concentrations (fugacities). 
The small (large) change of the local structure in g cc for large (small) polymers, respectively, is 
also seen in these simulations. Our results at higher polymer concentrations further show that 
the polymer induced depletion attraction is decreased by polymer-polymer excluded volume 
interactions as has been observed in neutron scattering experiments Jio| | . 

The structure of the colloidal liquid at triple coexistence has been the focus of a recent 
light scattering study which determined S cc (q) for 1.4 < qa c < 7.2 Q. Figure [| shows 
the data and theoretical curves evaluated without adjustable parameter. The small-q 1 colloid 
structure shows a strong dependence on polymer size. The theory accurately predicts the 
long wave length fluctuations which correspond to an order of magnitude enhancement of 
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Fig. 4 - Colloid structure factors compared to the experimental data by Moussaid et al. The 
parameters are (<j>c,(p P ,£c) = (0.333, 0.069, 0.085; short dashes and o), (0.404, 0.069, 0.131; long 
dashes and []), and (0.444, 0.053, 0.202; solid line and o). The inset shows the reduced polymer 
concentration dependence of the normalized colloid structure factors at the peak S p — S C c{q P )/ Qc, 
and at small wave vectors, So = S C c(0)/qc, at fixed colloid packing fraction, (f> c = 0.40, and for 
increasing polymer sizes, £ c = 0.03 (dot-dashed), 0.085 (short dashed), 0.202 (solid), and 1 (thin 
solid line). 




the osmotic compressibility relative to pure hard spheres. The large angle colloid structure 
experimentally shows little discernible polymer size dependence. The theory correctly captures 
the large-q trend concerning the peak intensity and location, except for the smallest polymer 
where a larger decrease in the peak height is predicted than seen in the experiments. The 
insets in fig. |J contrast the different trends of the colloid structure at small and large wave 
vectors when changing the polymer concentration for different polymer sizes. Whereas, the 
colloidal osmotic compressibility, S cc (q = 0)/g c , increases mainly due to approaching the 
spinodal, non-monotonic variations in the large angle scattering peak height, S cc (q p ), arise 
for small polymers corresponding to a polymer mediated distortion of the local collective 
packing ("cage") around a particle. Based on dynamic mode coupling ideas (2lj we expect 
this effect to cause the melting of the colloidal glass upon addition of small amounts of small 
polymer (£ c = 0.03) as reported in ||. 

In the limit where the spheres act as small depletants for the larger polymers, the descrip- 
tion of the intermolecular structure functions simplifies. Except for corrections of 0(a c /£ c ), 
the packing of the colloids becomes hard sphere like. The polymer segment-segment pair 
correlation function becomes long ranged, 

g pp (r) 1+ 2^( e ~^ -e"^) for £ c > a c , (6) 

where the collective correlation length, with £ c /£* = l+2<^ p (l— /oo(^c))) describes the decay 
of the polymer density from its (intermolecular) contact value, foo{4>c) = 2\i(i-4> ~)(i+2$) ' 
to zero. Polymer phase separation can be brought about by increasing the colloid density, 
as a spinodal instability (l/£* = 0) exists. This phase separation requires finite polymer 
concentrations, leads to a dense phase of interpenetrating polymers, and results from the 
increased intermolecular attraction induced by the spherical depletants. 
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To summarize, building upon the successful description of the structure of hard-spheres 
and athermal polymer solutions, we have developed a microscopic theory of the thermodynam- 
ics and structural correlations of binary particle-polymer mixtures. The inclusion of many- 
body interaction effects proves crucial for a successful description of the CPM-structures at 
higher densities. The explicit consideration of the conformational entropic contributions of 
the polymers is required to address the packing of larger polymers into the void space between 
colloidal spheres. Whereas solely entropic (steric) effects are considered extensions to include 
temperature dependent solvent quality effects and attractions have been achieved for dilute 
systems p2| , and can in principle also be incorporated for finite concentrations. 
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